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The influence of a wetting liquid on the coefficient of restitution (COR) is investigated experi- 
mentally by tracing freely falling particles bouncing on a wet surface. The dependence of the COR 
on the impact velocity and various properties of the particle and the wetting liquid is presented and 
discussed in terms of dimensionless numbers that characterize the interplay between inert ial, viscous, 
and surface forces. In the Reynolds number regime where the lubrication theory does not apply, 
the ratio of the film thickness to the particle size is found to be a crucial parameter determining the 
COR. 

PACS numbers: 45.70.-n, 45.50.Tn, 47.55. Kf 



I. INTRODUCTION 

The coefficient of restitution (COR), first introduced 
by Newton \^ as the ratio between the relative rebound 
and impact velocities of a binary impact, has been a sub- 
ject of continuous interest over centuries, along with the 
development of elastic Q , viscoelastic and plastic 
theories [H, Q. It characterizes the energy dissipation 
associated with the impact, which plays a key role in un- 
derstanding the collective behavior of macroscopic par- 
ticles, i.e. the dynamics of granular matter 0, This 
is largely due to the fact that the dissipative nature of 
granular matter arises from the inelastic collisions at the 
particle level. 

Due to its omnipresence in nature and various indus- 
tries, granular matter has drawn great attention from 
both physical and engineering communities in the past 
decades [9]. Concerning the modeling of granular matter, 
an appropriate collision model is essential for the success- 
ful implementation of kinetic or hydrodynamic theories 
to granular matter [10|-[l3|, see for example the dynam- 
ics of Saturn's rings [1^, or the pattern formation under 
vertical agitation [15|. Despite those successful exam- 
ples for dry granular matter, a continuum description for 
wet granular matter, which considers the cohesion arising 
from the wetting liquid phase, is still far from established 
[HI, [13 • Therefore, in order to provide a solid basis for 
a continuum modeling of wet granular fiow - for exam- 
ple to describe natural disasters such as debris fiow - a 
thorough understanding of the dynamics associated with 
wet impacts is desirable. 

With the development of pharmaceutics, mining and 
food industries, the COR for wet impacts has become an 
important issue for the engineering community in terms 
of decoding the underlying physics associated with the 
agglomeration of particles with liquid binders. The pi- 



oneering work by Rumpf [l8| half a century ago has in- 
cluded a detailed description of the capillary force of a 
pendular bridge and treated it as the dominating cohesive 
force in determining the statical behavior of wet granular 
matter, e.g. the tensile strenth [l9|. Later on, the viscous 
force has been found to play an important role in typical 
granulation processes, too |2Q|-[27j. And a dynamic liquid 
bridge could be an order of magnitude stronger than a 
quasi-static one [28|, [29| . Binary impacts of particles with 
viscous liquid coating have been extensively investigated 
by experiments and models using lubrication theory (3Q| - 

Recently, a new class of collision experiments in- 
vestigating three-body interactions utilizing a "Stokes's 
Cradle" apparatus, revealed a more complex scenario 
(35I , [36| : In addition to the agglomerated and separated 
final states known from binary impacts, partially sepa- 
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FIG. 1. (color online) Sketch of the experimental setup. The 
bouncing of the glass bead, initially held by the vacuum noz- 
zle, on the glass container is recorded with a high speed cam- 
era (Camera 2). The thickness of the liquid layer is monitored 
by detecting the laser beam (red line) reflected from the liquid 
surface with Camera 1. 
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FIG. 2. A series of snapshots captured with a frame rate of 450 Hz showing a 4 mm glass bead bouncing on ; 
with a 1 mm water film. 



glass plate covered 



rated outcomes including a "Newton's cradle" and a "re- 
verse Newton's cradle" collision could be observed. 

Despite all those investigations, a well tested collision 
law suitable for modeling the dynamics of wet granular 
behavior (STMsof , as well as a comprehensive knowledge of 
the energy dissipation associated with the impact, is still 
lacking. In the current work, the COR of a ball bounc- 
ing back from a flat lubricated surface is investigated as 
a function of the impact velocity and various particle as 
well as liquid properties. From this, the kinetic energy 
dissipated during the impact process is derived and dis- 
cussed within the framework of existing models. 



II. EXPERIMENTAL SETUP AND 
PROCEDURE 

Figure [1] shows a sketch of the experimental setup 
used for the COR measurements. Spherical glass beads 
(SiLiBeads type P) with a diameter range from D = 
2.8 mm to 10 mm, roughness ~ 5 /im, and density pg = 
2.58 g/cm^ are used in the experiments. By controlling 
the pressure in the vacuum nozzle, the initially wet par- 
ticle is released from its starting position and falls onto 
a glass container (20 cmx 5 cm) covered with a thin layer 
of wetting hquid: Purified water or silicone oil (Carl 
Roth M5 and M50) with dynamic viscosity rj of 4.6 mPas 
and 48mPas, respectively. The bottom of the container 
is 2 cm thick. It is leveled within 0.03 degrees, so that 
bouncing on various positions in the container explores 
a similar liquid layer thickness. The layer thickness 5 is 
measured by detecting the shift of a laser beam reflected 
from the surface of the liquid and the glass plate with a 
CCD camera (Camera 1, Lumenera Lul35). The mirror 
attached to the bottom of the container creates multiple 
reflections of the laser beam, in order to enhance the sen- 
sitivity of the device. The length of the mirror (7.8 cm) 
is chosen as a compromise between the sensitivity and 



the field of view. By fixing the container, laser and the 
camera on a leveled optical table, the error of the film 
thickness measurement could be minimized to a satisfac- 
tory level (< 10 /im). 

To obtain the impact and rebound velocities, the 
bouncing of the particle is recorded by a fast camera 
(Photron Fastcam Super lOK) with succeeding digital 
image processing of individual frames. A close view of 
the colliding event, as shown in Fig. [21 clearly demon- 
strates the important role that the wetting liquid plays 
during the impact. As the sphere hits the liquid sur- 
face, a circular wave front occasionally accompanied with 
a splash will be generated. As the ball rebounds from 
the surface, a liquid bridge will form between the sphere 
and the hquid surface, which continuously deforms and 
elongates until it ruptures at a distance larger than the 
particle diameter. Associated with the rupture event, 
satellite droplets may form, which bounce on the liq- 
uid surface and coalesce partially into smaller droplets 
[iol , IZlj . Obviously, the formation of wave fronts, defor- 
mation and rupture of liquid bridges, the viscous force, 
and the added mass to the sphere due to the wetting hq- 
uid will all contribute to the mechanical energy reduction 
of the impacting particle, which in turn leads to a smaller 
COR compared with dry impacts. 

Figure [3] illustrates the influence of wetting by pro- 
viding a comparison between the trajectories obtained 
from wet and dry impacts. The particle diameter is 
D = 5.5 mm, and the film thickness of the silicone oil 
M5 is (5 = 225 /im in the wet case. To determine the 
center of the sphere, the image processing procedure em- 
ploys a Hough transformation [42| (upper panel of Fig. [3]). 
Subsequently, each bouncing trajectory is extracted and 
subjected to a parabolic fit (see the solid line in Fig.[3l^a) 
as an example), in order to obtain the peak position /ipeak 
and the impact velocity. If the normal COR (en) is inde- 
pendent on the impact velocity, the velocity after the zth 
rebound will be related to the first impact velocity '^o by 
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FIG. 3. (color online) Trajectories of a particle bouncing on 
a wet (a) and a dry (b) surface after image processing. The 
image sequence in the upper panel represents a fraction of 
the wet trajectory (a) with superimposed centers and bound- 
aries of the sphere. The blue line in (a) corresponds to a 
parabolic fit to the trajectory after the first rebound. The 
peak positions of the trajectories /ipeak obtained from the fits 
are marked with triangles in (a) and upside down triangles in 
(b). The inset in (a) shows /ipeak as a function of the number 
of impacts. 



'^i = e^^o- This leads to a linear decay of the peak height 
/ipeak with the number of impacts z in a semi- log plot, 
according to Ig/ipeak = Ig/^o + 2zlgen, with hpeak oc vf. 
The initial falling height ho and Cn determine the offset 
and slope of this line. As shown in the inset of Fig.[3l^a), 
the logarithm of /ipeak decreases linearly with the number 
of impacts for dry impacts, indicating that the normal 
COR stays almost constant for the number of impacts 
measured here. In a recent work on dry impacts [43|, 
a more detailed analysis reveals that the dry COR de- 
creases slightly with the increase of ^'impact- However, 
this dependence is much weaker than the one for wet im- 
pacts, on which we are focusing here. In this case, the 
variation of the slope indicates that the COR for wet im- 
pacts decreases strongly with the number of impacts, i.e. 
with the impact velocity. 



III. EXPERIMENTAL RESULTS 

As shown in Fig. HI this dependence of the COR on 
the impact velocity '^impact is qualitatively the same for 
various particle diameters and for both silicone oil (M5) 
and water wetting: It grows initially with 'i^impact and 
saturates at a certain value varying from 0.8 to 0.9. For 
a fixed '^impact, the COR decreases systematically with 
particle diameter for both silicone oil and water wetting. 
The monotonic growth of the COR on D reveals that the 
influence of the wetting liquid decreases as the inertia of 
the particle (oc D^) increases. In the lower panel of Fig. 21 
the rebound velocity '^^rebound is plotted as a function the 
impact velocity 'i^impact- Similar to the case without wet- 
ting liquid (shown as a gray dashed line), 'i^rebound grows 
linearly with ^impact, and the trend is qualitatively the 
same for all parameters used here. Therefore, fitting the 
data with '^rebound = einf(vimpact - ^c) givcs risc to two 
parameters that characterize the impact velocity depen- 
dence: A slope einf corresponding to the COR at infinite 
'^impact, i-e. the saturated value of the COR, and an off- 
set Vc corresponding to a critical energy Ec below which 
no rebound would occur. Ec = mv'^/2 is obtained from 
the intersection Vc of the linear fits shown in the lower 
panels of Fig. 2] with the x-axis, where m is the mass of 
the particle. 

Figure [5] shows the corresponding results. For the dry 
impacts, the results stem from linear fits of the data ob- 
tained for different particle sizes. The critical energy for 
the dry impacts is very close to 0. In contrast, the crit- 
ical energy for wet impacts is on the order of few /iJ. It 
shows a monotonic decay for water wetting, and a more 
complicated relationship for M5. As shown in Fig.[5](b), 
einf - the upper limit of the COR - varies from 0.8 to 
0.9, and is generally smaller than edry From the energy 
dissipation perspective, this demonstrates that the ratio 
between the energy dissipation from the wetting liquid, 
A^wet , and the kinetic energy at impact £^impact will not 
diminish as '^impact grows. For both silicone oil (M5) and 
water wetting, einf shows similar values with weak de- 
pendence on the particle sizes, although M5 silicone oil 
is 5 times more viscous than water. 

The above experimental results indicate that the COR 
depends strongly on the impact velocity as well as var- 
ious particle and liquid properties. In order to explore 
the relation between COR and all these parameters, it is 
essential to have a proper classification of the parameters 
in terms of dimensionless quantities that characterize the 
relation between inertia, viscous and capillary effects. In 
the case where the viscous force dominates, the dynamics 
of wet impacts could be described by the elastohydrody- 
namic theory [2l|, [sO, [H, [s^. Within this framework, 
the Stokes number, defined as St = pgDvimpact/^V^ is 
usually used to characterize the ratio between the inertia 
of the particle and the viscous force. Normally, this case 
is justified by the criterion that the Reynolds number of 
the wetting liquid Re <^ 1 [s^, where Re = '^impact Z^, 
with pi the density and 5 the thickness of the wetting liq- 
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FIG. 4. (color online) Normal restitution coefficient Cn and rebound velocity t^rebound as a function of impact velocity t;impact 
for the impacts of particles with various diameters D on silicone oil (left column) and water (right column) wetting surfaces 
with fixed film thickness 6=1 mm. The solid lines in the lower panels are linear fits to the data and their representatives are 
shown in the upper panels as a guide to the eyes. The dashed gray lines in the upper panels represent the normal restitution 
coefficient Cdry = 0.976 db 0.001 for dry impacts, which is obtained by a linear fit of the data for all particle sizes (gray diamonds 
shown in b). Error bars smaller than the symbol size are not shown. 



uid. This implies the conditions that either the liquid is 
highly viscous, or the film thickness is small enough for 
the lubrication theory to apply. Within this limit, the 
contribution from the wetting liquid to the total energy 
dissipation is mainly due to viscous damping. Although 
the range of Reynolds numbers for the current investiga- 
tion (10 < Re < 10^) suggests that the role that the vis- 
cous force plays may not be prominent, we represent the 
COR results shown in Fig. |4] with respect to the Stokes 
number as a starting point. 

As shown in Fig.[6l the scaling of with 'i^impact and 
D for more viscous silicone oil (M5) could be better de- 
scribed by the Stokes number St than that for water. 
For M5 wetting (corresponding to Re = 20-360), data 
for various D show a general trend of initial growth from 
St ^ 100 to 500, followed by a saturation to einf between 
0.8 and 0.85. Concerning the case of water wetting (cor- 
responding to Re = 100-1800), the scatter of the data 
obtained with various particle sizes (shown in Fig.[6l^b)) 
is much more prominent than for the case of M5 wetting. 
Although the trend of a significant growth followed by a 
saturated value persists, both the slope of increase and 



the saturated value differ as D varies. This difference be- 
tween wetting liquids with various viscosity suggests that 
further parameters have to be considered as the Reynolds 
number goes beyond the regime where the elastohydro- 
dynamic theory applies. 

From another point of view. Fig. [6] also reveals a rela- 
tively small difference of the COR between M5 and water 
wetting, even though the corresponding viscosity ratio is 
5. This result clearly demonstrates that the COR is also 
determined by other liquid properties. As an example, 
the surface tension of water (7 = 72.0 mN/m) is much 
larger than that of M5 (19.2 mN/m), which may lead to 
a larger energy dissipation from the formation of capillary 
waves and the break of capillary bridges upon rebound. 
To focus on the effect from viscosity, we use wetting liq- 
uids with similar surface tensions for comparison in the 
following part of the paper. 

The influence of the liquid film thickness 6 and the 
dynamic viscosity r] on the wet impacts is presented in 
Fig. [71 The upper part (a) shows the corresponding nor- 
mal restitution coefficient en = '^rebound /'^impact as a func- 
tion of '^impact- For wet impacts, shows clearly a de- 
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FIG. 5. The critical energy Ec and the saturated value of the 
COR Cinf as a function of particle diameter D. The solid line 
at = is a guide to the eyes. 
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pendence on '^impact, and this dependence is more promi- 
nent for a liquid film with higher viscosity. With M5 as 
wetting liquid, Cn grows slightly with 'Z^impact and tends 
to saturate at a value depending on the film thickness 
S. Similar to the results shown in Fig. HI this value is 
smaller than edry, revealing the extra energy dissipation 
from the wetting liquid. As the film thickness 5 increases, 
Cn decreases systematically for both wetting liquids, be- 
cause the viscous damping force is effective over a larger 
distance. A comparison between the influence from the 
viscosity and the film thickness indicates that the former 
has a stronger impact on the energy dissipation and en- 

In Fig.[7Jb), the relation between '^rebound and '^impact 
is presented. For dry impacts, a linear growth of the re- 
bound velocity 'i^rebound with the impact velocity 'i^impact 
without an offset suggests a constant edry = 0.985 for the 
range of 'i^impact explored here. For wet impacts, 'i^rebound 
decreases systematically with the liquid film thickness S 
at a certain '^impact- As the liquid viscosity increases by 
an order of magnitude (from M5 to M50), this trend is 
more prominent, indicating the crucial role played by the 
viscous damping. The growth of 'i^rebound with 'i^impact 
could well be fitted with a straight line. The fit again 
gives rise to a slope einf that is smaller than edry and a 
threshold energy Ec below which no rebound would oc- 
cur. As shown in Fig. [H^a), this threshold is, for M50, 



FIG. 6. (color online) Normal restitution coefficient Cn as a 
function of the Stokes number St for both silicone oil (a) and 
water (b) films. Parameters are the same as in Fig. |4l 



more than an order of magnitude larger than that for M5, 
within the range of film thicknesses used. This suggests 
that Ec depends not only on the rupture energy of a cap- 
illary bridge, but also on viscous damping. As shown in 
Fig.[8fb), the slope einf is not strongly influenced by vis- 
cosity compared with Ec. For relatively thin film wetting, 
en could be the same within the error bars. The slope 
einf stays constant within the range of film thickness and 
decays shghtly for more viscous M50 wetting. 

Since the influence of the viscosity could be charac- 
terized by the Stokes number, we replot en, now as a 
function of St. As shown in Fig.[9l^a), the COR grows 
dramatically at small St^ which corresponds to the data 
of the M50 wetting case, and saturates at larger St. For 
various film thickness 6, this trend is qualitatively the 
same. Quantitatively, the saturated value einf decreases 
as S grows, indicating that the assumption of small di- 
mensionless film thickness 6/D for the lubrication theory 
to apply could not be fulfilled for the whole range of St in- 
vestigated here. In other words, this feature shows that a 
further dimensionless parameter including the film thick- 
ness should be considered beyond the viscous dominating 
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FIG. 7. (color online) The normal restitution coefficient Cn 
(a) and rebound velocity t^rebound (b) as a function of impact 
velocity for impacts of a glass bead with D — 5.5 mm 

on dry and wet surfaces covered with silicone oil M5 and M50. 
S denotes the film thickness. The error bars correspond to the 
statistical error over 10 runs of experiments for the wet im- 
pacts. Solid lines in (b) are linear fits to the corresponding 
data. Their representatives are shown in (a) as a guide to 
the eye. Squares and dots correspond to M5 and M50, re- 
spectively. For dry impacts, the restitution coefficient Cdry is 
0.985 with an error of 0.001. 



regime. 

This parameter is chosen as the dimensionless length 
scale S = S/D, because it ties the Stokes number with 
the Reynolds number of the wetting liquid. According to 
this definition, the ratio between the Reynolds number 
and the Stokes number is Re/ St = 96 where p = p\/ pg 
is the density ratio between the liquid and the particle. 

In Fig.[9l^b), further experiments with the restriction 
6 ~ 0.04 are presented. In contrast to Fig.[9ja), the data 
from various film thicknesses coincide over a wide range 
of St if 6 is fixed. It also gives rise to a master curve 
= einf(l — Stc/St)^ as indicated clearly in the inset. 
The linear fit yields einf = 0.908 ± 0.002, and a critical 
Stokes number Stc = 14.00 ± 0.20. Therefore, the usage 
of the Stokes number as a control parameter could be 
extended to the regime Re > 1 and large film thickness, 
provided that the dimensionless length scale 6 is kept 
constant. 
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FIG. 8. The critical energy Ec and the saturated value of the 
COR einf as a function of the film thickness S. Ec is obtained 
from the intersection Vc of the linear fits shown in the lower 
panel of Fig. [71 with the x-axis. Cinf corresponds to the slope 
of these fits. 



IV. ANALYSIS OF THE ENERGY DISSIPATION 

To understand the dependence of the COR on various 
particle as well as liquid properties, it is helpful to ana- 
lyze the associated energy dissipation. With £^diss being 
the total kinetic energy loss of the particle during the im- 
pact, the dependence of the COR on the kinetic energy 
at impact Ei can be written as 



\/l — Ediss/Ei. 



(1) 



The dissipated energy £^diss can be treated as the sum 
of two parts: the part transferred into the solid body 
A£^dry, and the other part taken by the wetting liquid 
A£;wet, i.e. 



dry 



•AK 



(2) 



Provided that the two parts are independent from each 
other, i.e. the wetting liquid does not change the energy 
dissipation of the dry impact, AE^^t could be obtained 
experimentally by 



dry 



el). 



(3) 
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FIG. 9. (color online) The normal restitution coefficient as 
a function of the Stokes number St. (a) corresponds to the 
data shown in Fig|7l (b) corresponds to the data with the 
dimensionless film thickness 6/D roughly constant. The error 
bars correspond to the statistical error of 10 runs of experi- 
ments. The solid curve in (b) corresponds to the master curve 
Gn = 0.908(1 — 14.00/5't), which is obtained from a linear fit 
to all the data shown in the inset. 



The whole process of the colliding event can be sepa- 
rated into two parts: impact and rebound. During the 
impact, the kinetic energy of the particle will partly be 
transferred to the wetting liquid. This amount of energy 
will finally be dissipated by the motion of the viscous 
liquid, including surface waves or even splashes, depend- 
ing on the competition between the inertial, viscous and 
surface forces. During the rebound, the rupture of the 
capillary bridge will lead to a certain amount of surface 
energy loss in addition to the damping caused by the 
motion of the liquid. Moreover, the mass of the wetting 
liquid dragged away by the sphere might lead to a further 
reduction of the COR. Based on the above analysis, one 
can take the most prominent terms and use 



(4) 



to estimate AE'wet theoretically , where A£^visc represents 
the energy dissipated via the viscous damping force act- 
ing on the particle, AE'b corresponds to the energy loss 
arising from the surface energy change of the fluid, and 
A^acc is the kinetic energy change of the fluid before and 



after the colliding event. 

In the limit that the thin fllm lubrication theory ap- 
plies, the viscous force acting on the particle can be es- 
timated by = 37rr]D'^Vimpact/'^x [33|, where x denotes 
the distance between the sphere and the plate. Follow- 
ing Ref. [2l|, one might assume the same force law for 
both approach and departure of the sphere. By integrat- 
ing over the distance that the viscous force applies, we 
obtain 



3 S S 

AE^isc = 2^^^^ ^impact (In - + In -^), (5) 

where e = 5 /im is the roughness of the sphere, and 
is the rupture distance of the liquid bridge. For a crude 
estimation, we take a flxed = 21) according to the 
snapshots taken and assume that the velocity does not 
change during the impact. 

In Fig.[TOl AE^et for the experimental results shown 
in Fig. [71 is plotted in comparison with A£^visc- Quali- 
tatively, the monotonic growth of the energy dissipation 
with the impact velocity, and the increase of energy dis- 
sipation with the fllm thickness agree with the estima- 
tion from Eq. ([5]). However, quantitatively a comparison 
between the estimated viscous damping term A£^visc and 
AE^Qt reveals that a substantial amount of the latter can 
be attributed to the viscous damping for the case of M50 
wetting, while this term plays a much weaker role for the 
case of less viscous silicone oil M5 wetting. Obviously, 
the growth of the difference with increasing 'i^impact and 
decreasing viscosity represents a deviation from the limit 
that lubrication theory applies. Taking 30% deviation as 
the limit, one can estimate the corresponding Reynolds 
number to be Re ~ 10 for the case of 5 « 0.04. For M5 
wetting, the growth of the energy dissipation A^^et with 
'^impact deviates slightly from a straight line, suggesting 
that the dominating energy dissipation term has a higher 
order dependence on the impact velocity. 

The second term in Eq. ^ stems from capillary forces. 
Upon rebound of the sphere, a liquid bridge may form 
between the sphere and the liquid surface. The corre- 
sponding energy dissipation due to the deformation and 
rupture of this liquid bridge can be estimated by an in- 
tegration of the force arising from the surface tension 
over the length that it acts. This capillary force has two 
components: the surface tension acting on the perime- 
ter of the neck (27rrn7 with Vn the neck radius and 7 
the surface tension), and the second part arising from 
the Laplace pressure pb that acts on the cross-section of 
the neck (— ApbTrr^). Based on quasi-static experimental 
veriflcations, a close form approximation of the capillary 
force Fc between two spheres has been given as 



7TDjC0s{(j)) 

1 + 2.16'* + 105*2' 



(6) 



where 6* = s^^ D/2V]j is the half separating distance s 
rescaled by the characteristic length scale ^JD|2V\^ with 
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FIG. 10. (color online) Energy dissipation A^wet due to wet- 
ting as a function of impact velocity with the wetting liquid 
properties the same as in Fig|7]for both silicone oil M50 (a) 
and M5 (b) wetting. Solid lines represent estimated values 
of the energy dissipation A^^visc from viscosity (see text for 
detailed descriptions). The line colors are in accordance with 
the data points for various film thickness. 

the bridge volume Vb, and (j) corresponds to the contact 
angle |44| . 

Taking the rupture distance 5^ as the integration hmit, 
one could estimate the rupture energy of the liquid bridge 
to be 



A^b « T^iV^m^. (7) 

A rough estimation of the bridge volume Vb ^ 1^^/16, 
based on the snapshot taken, gives rise to ^ 0.7/iJ for 
silicone oil wetting a glass bead with diameter 5.5 mm. 
Considering the energy dissipation obtained by the COR 
measurements shown in Fig.[10l plays a minor role for 
the few mm sized particle used here. Note that plays 
a more prominent role as D decreases, because its growth 
with is in contrast to £^visc D^. 

As demonstrated in Fig(TOl both the damping from 
the viscous force and the rupture of liquid bridges can 
not explain the amount of energy dissipation for the case 
of silicone oil M5 wetting. Therefore, other effects, like 
e.g. the inertia of the liquid or surface waves, should be 
considered. 

As a first approximation, the inertial effect could be 



FIG. 11. (color online) Data points are the same as shown 
in Fig. 1101 Solid lines represent estimated values with the 
consideration of both viscous damping A^^visc and the energy 
transfer to the fluid AEacc (see text for detailed descriptions) . 
The line colors are in accordance with the data points for 
various film thickness. 

estimated from the kinetic energy of the liquid being 
pushed aside by the impact [45|- The volume of the 
hquid can be estimated by the spherical cap immersed 
in the liquid film V = 7rD^P{l/2 - S/3). From the 
length scale taken as the base radius of the spherical 

cap -^/l — (1 — 2S)'^D/2 and the time scale ^/'i^impact for 
the particle to penetrate the liquid layer, one estimates 

the average velocity vi = '^impact 1/S — 1. As a con- 
sequence, the kinetic energy A£^acc of the liquid being 
pushed aside yields 

A^acc = IpiVvf = 3p{S - + p^)Eu (8) 

which shows a linear dependence on the kinetic energy Ei 
of the impact particle. Figure [TT] shows that, by taking 
both A£^visc and A£^acc into account, the influence from 
the inertia effect is more prominent for less viscous M5 
wetting. The combination of both forces leads to a better 
agreement with the experimental data, when compared 
to Fig. fTOT b). However, considering both the inertial and 
the viscous damping parts of the energy dissipation can 
not explain the experimental results for less viscous M5 
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wetting quantitatively. This indicates that further the- 
oretical considerations, e.g. on additional energy dissi- 
pation terms, or a more careful characterization of the 
inert ial effects, are desirable. 

The fact that the ratio between A£^acc and is not 
velocity dependent suggests that the inertia of the wet- 
ting liquid will not contribute to the impact velocity de- 
pendence of the en- It does, however, explain why einf 
obtained from linear fits of the data is generally smaller 
than edry Based on the Eqs. ([3j) and (j4|), the Ei depen- 
dent COR could be written as 



In the limit of large '^^impact, the last term AE\)/Ei 
could be safely ignored so that two independent param- 
eters are enough to determine the impact velocity de- 
pendence of the COR. A comparison to the linear fits 
en = einf('i^impact " '^c) employed before immediately re- 
veals that the linear fit is a first order approximation of 

the Eq. (j9j) and einf = ^J^^^^Xe^^^Je-,. By ignoring 

the higher order terms of 5 in A£^acc7 one derives a lin- 
earized form einf = ^dvy — 3p(5/2edry, which suggests a 
reduction of einf from edry by about 0.06 for typical ex- 
perimental values of p = 1/2.5 and ^ = 0.1. A compari- 
son to Fig.lSl^b) shows that this simplified formula gives 
the correct order of magnitude for the einf estimation. 
Moreover, the monotonic decrease with 5 is captured by 
this formula qualitatively, except for the 75 fiia thick M5 
wetting case. 



V. CONCLUSION 

In summary, the normal coefficient of restitution 
(COR) for a free falling sphere on a wet surface is in- 
vestigated experimentally. The dependence of the COR 
on the impact velocity and various particle and liquid 
film properties is discussed in relation to the energy dis- 
sipation associated with the impact process. 

i) For dry impact, the COR corresponds to the slope 
of the rebound vs. impact velocity. For wet impacts, the 
rebound velocity and the impact velocity are also found 
to fall into a straight line, but with a smaller slope and 
an offset corresponding to a finite critical impact velocity. 
Even though linear fitting is only a first order approxima- 
tion of en, it successfully characterizes the impact veloc- 
ity dependence of the COR with two parameters einf and 
E^. Therefore this simplification is justified to be a good 
candidate for computer simulations aiming at modeling 
wet granular dynamics on a large scale. 

ii) The dependence of the COR on the impact veloc- 
ity, dimension of the sphere and the viscosity of the liquid 
could be well characterized by the Stokes number, which 
is defined as the ratio between the inertia of the sphere 
and the viscosity of the liquid, provided that the dimen- 
sionless length scale 5 is fixed. This result supports the 
usage of the Stokes number for scaling the data, even 
beyond the low Reynolds number regime where it has 
originally been introduced. 

iii) Concerning the energy dissipation arising from the 
wetting liquid, the viscous damping term dominates for 
Reynolds number up to i?e ^ 10. Away from that limit, 
further effects, such as the inertia of the liquid film, have 
to be considered. The rupture energy of a capillary bridge 
during the rebound process could be safely ignored for 
the few mm sized particles used here. 
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Ennis and Jiirgen Vollmer for helpful hints. This work 
is partly supported by the German Science Foundation 
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plexer Kontinua' through Grant No. Kr 1877/3-1. 
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